The phase transition of triplet reaction-diffusion models 
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The phase transitions classes of reaction-diffusion systems with multi-particle reactions is an 
open challenging problem. Large scale simulations are applied for the 3A — > 4A, 3A — > 2A and the 
3A — » AA, 3A —> triplet reaction models with site occupation restriction in one dimension. Static 
and dynamic mean-field scaling is observed with signs of logarithmic corrections suggesting d c — 1 
upper critical dimension for this family of models. 
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I. INTRODUCTION 

The exploration of the universal behavior of non- 
equilibrium phase transitions is still an open problem 
of statistical physics 0, 0, 0- In equilibrium critical 
phenomena symmetries, dimensions and boundary con- 
dition are the relevant factors determining a universal- 
ity class. During the study of non-equilibrium models 
other circumstances, like initial conditions or topological 
effects in low dimensions have been proven to be deci- 
sive . Unfortunately solid field theoretical background 
exists for only a few reaction-diffusion (RD) systems ex- 
hibiting continuous phase transition to absorbing states 
|{|. These are mainly branching and annihilating ran- 
dom walk models A — > (n + I)A, 2A — > built up from 
unary particle creation reactions 0, 0, @ RD models are 
interesting since many other types of systems like surface 
growth, spin systems or stochastic cellular automata can 
be mapped on them [3| . 

Recently RD systems with multi-particle creation has 
became in the focus of research. Numerical studies re- 
sulted in debated critical phenomena and generated a 
long series of publications. An interesting example was 
being investigated during the past decade that emerges 
at phase transitions of binary production systems [l| 
(PCPD). In these systems particle production competes 
with pair annihilation and single particle diffusion. If the 
production wins steady states with finite particle density 
appear in (site restricted) models with hard-core repul- 
sion, while in unrestricted (bosonic) models the density 
diverges. If the annihilation is stronger an absorbing 
phase emerges, which is either completely empty or con- 
tains a solitary diffusing particle. In between the two 
phases a continuous phase transition can be observed in 
site restricted models. 

In triplet reaction models at least three particle is 
needed to contact for a reaction. They have been in- 
vestigated by simulations E(l ITTl IT^ and numerical 
Langevin equation solution 15J . The first simulation re- 
sults [hJ [13 for 3A -»■ 4A, 3A -» 2A models in one 
dimension claimed a distinct universal behavior from the 
known ones. Simple power-counting analysis of an ef- 
fective Langevin equation corresponding to the coarse 
grained microscopic model results in d c = 4/3 How- 



ever the numerical estimates for the critical exponent 
describing the density decay from homogeneous random 
initial state p(t) cx t~ a differed significantly a — 0.32(1) 
[Tof vs a = 0.27(1) In the former case a site re- 

stricted model was investigated and scaling was reported 
in the 10 4 < t < 10 6 region. In different suppressed 
bosonic triplet models - where the multiple site occu- 
pancy is suppressed by an exponential probability factor 
- exhibited scaling for 10 4 < t < 10 7 Monte Carlo steps 
(MCS) (throughout the paper the time is measured by 
MCS). Renormalization group analysis pointed out 
that a single field theory does not exhibit a nontrivial 
stable fixed point and suggested d c = 1 for such models. 
This study raised the possibility that a proper field the- 
ory should be a coupled one, with positively correlated 
clustered particles and solitary random walkers. An other 
simulation study [l^ on site restricted models reported 
scaling agreeing with mean-field exponents a — 0.95(5) 
and (3 = 1.07(10), where (3 is the order parameter expo- 
nent in the active phase p cx \p — p c \^ ■ 

Very recently a coupled field theoretical description of 
such systems is suggested |l5|. An effective Langevin 
equation between a DP like and an annihilating random 
walk (ARW) system is analyzed by numerical integra- 
tion technique. Note however the ARW like system is 
described by a positive noise term saying that at the 
critical point the anti-correlations do not play a role and 
serve merely as a fluctuating source to the primary field. 
As the consequence this field theory leads to the same 
critical scaling behavior as that of the PCPD albeit with 
a different upper critical dimension d c = 4/3 (vs d c = 2 
for PCPD). So according to this study in one dimension 
one should see PCPD class exponents: a — 0.20(1) and 

^ = 0.40(1) mill. 

In the present study I extend the simulation time of the 
3A — » 4A, 3A — > 2A model investigated in [lj] by two 
orders of magnitude and follow the static and dynamic 
scaling behavior of particles and triplets (AAA) at differ- 
ent diffusion probabilities. In Section (jlllll I apply this 
the same kind of analysis for the 3A — > 4A, 3A — > 
model. 
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FIG. 1: (Color online) Density decay of the 3A -> 4A, 3A -> 
2A model at D — 0.1. Different curves correspond to p — 
0.301, 0.302, 0.3027,0.3031, 0.3032, 0.3033, 0.30335, 0.3034, 
and 0.3035 (from top to bottom). The insert shows the decay 
of triplets: p 3 t 1/3 , for p = 0.30337 



FIG. 2: (Color online) Density decay of the 3A -> 4A, 3A ~* 
2A model at D = 0.8. Different curves correspond to p = 
0.406, 0.407, 0.4074, 0.4075 0.408, 0.409 (from top to bottom). 
The insert shows the triplet density decay at p = 0.4075. 



II. SIMULATIONS OF THE 3,4 
MODEL 



4A, 3A -> 2A 



The simulations were carried out on 2 x L = 10 5 sized 
systems with periodic boundary conditions. The initial 
states were randomly half filled lattices, and the density 
of particles, singlets and triplets is followed up to 2 x 10 9 
MCS by random sequential dynamics. An elementary 
MCS consists of the following steps. A particle A is cho- 
sen randomly and the following processes are done: 

(a) A% <-» A with probability D, 

(b) 3 A -> 2 A with probability p(l - D), 

(c) 3A -> AA with probability (1 - p)(l - D), 

such that the reactions were allowed on the left or right 
side of the selected particle strings randomly. The time 
is updated by p[t). The time - measured by MCS — is 
updated by l/np, where np is the total particle number 
at time t. In order to get precise p(t), critical point and 
exponent estimates the number of independent realiza- 
tions varied between 20 and 350 per p and D throughout 
this study. 

First I extended the simulations at D = 0.1 published 
in [13 from t max = 10 7 MCS by a factor of 200 in 
time. Figure ^ shows the density decay multiplied by 
i 1 / 3 . Following an long initial transient, where the de- 
cay is slow an excellent agreement with the mean-field 
scaling can be observed for 10 6 < t < 10 9 MCS for 
p c = 0.30337(2). The triplet density decays in the same 
way (p3 oc t 1 / 3 ) as the total density suggesting d c = 1. 

I repeated the simulations for a higher diffusion rate 
D = 0.8. In this case due to the dynamics of this model 
the reaction rates are smaller and one can observe a faster 
than mean-field decay for 10 4 < t < 5 x 10 8 MCS. Agree- 
ment with mean- field scaling sets in for 5 x 10 8 < t < 



2 x 10 9 MCS steps for p = 0.4075(1). Fast transient de- 
cay for intermediate times has already been explained in 
case of the AA — > A coagulation model with finite re- 
action rates 0, Q] . One may expect similar behavior 
for the triplet annihilation case where the low reaction 
rates at D = 0.8 can explain the fast transient seen here. 
Again the triplet density decays by the mean-field law for 
t > 5 x 10 s MCS as expected at the upper critical dimen- 
sion d c = 1. Using the critical point estimates obtained 
by the dynamical simulations I investigated the singular 
behavior of the order parameter in the supercritical re- 
gion as well. The order parameter is expected to scale 
as 
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To get the steady state densities I followed the decay in 
several samples until saturation occurs and averaged it in 
a long time window exceeding a level-off seen on log-lin. 
scale. The data are analyzed by the local slopes method. 
According to this 
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and as pi — > p c the effective exponent tends to the true 
critical value j3 e ft — > (3. As Figure |21 shows the j3 e ff con- 
verges to the mean-field value j3 = 1 for both diffusion 
rates I investigated. For D = 0.8 on can see an over- 
shooting, while for D = 0.1 local slopes approach the 
asymptotic value from below following an upward cur- 
vature. This kind of effective exponent behavior usually 
corresponds to logarithmic correction to scaling. In case 
of the general form 



p(oo) = e^ln^e) , 



(3) 
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FIG. 3: (Color online) Effective static order parameter ex- 
ponent results in the active phase. Bullets correspond to 
D = 0.8, boxes to D = 0.1 of the 3 A -» 4A, 3 A -» 2 A 
model. Stars denote the results of the 3A — > 4A, 3A — > 
model. 



FIG. 4: (Color online) Density decay in the 3A -* 4A, 3A -> 
model at D = 0.8. Different curves correspond to p = 
0.1185, 0.1187, 0.11875,0.1188, 0.11886, 0.1189, 0.119 (from 
top to bottom). The insert shows the decay of triplets for 
p = 0.11886. 



where e = \p — p c the effective exponent behaves as 

dln(p(oo)) a; 
A " = dln( C ) =/?+ InW- (4) 

Applying this form for the data of Fig|S] one obtains: 
(3 = 1.05(5), x = 1.2(1) for L> = 0.1 and (3 = 0.93(10), 
x = —1.1(1) for D = 0.8. Although the assumed scaling 
correction form may look somewhat ad hoc it indicates 

4 = 1- 

I also considered the scenario suggested by fi^| accord- 
ing to which PCPD scaling (a = 0.19(1) HEEj) should 
be observed at the transition point. By assuming that 
such density decay appears for very long times (t > 10 9 
MCS) one can read-off the corresponding slightly differ- 
ent critical point estimates: p c = 0.30325(5) for d = 0.1 
and p c — 0.4071(1) for D = 0.8. Using these values in 
the local slopes analysis one can obtain: (3 = 0.85(5) for 
D = 0.1 and (3 = 0.80(5), neither of them is near to the 
exponent (3 = 0.40(1) of the PCPD class [13 • Therefore 
these simulations can't support the PCPD class scenario. 



III. SIMULATIONS OF THE 3A -> 4A, 3A -> 
MODEL 

Considering the mean-field type of scaling behavior 
of the 3A — ► 4A, 3A — > 2A model one may speculate 
that in this model the spatial fluctuations are somewhat 
suppressed: it's hard to create nearly arbitrarily large 
regions void of particles from a place where "annihila- 
tion" reactions have taken place due to the 3A — > 2A 
rule. In the low-diffusion regime, when a particle was cre- 
ated from the configuration ...0AAA0... it is most likely 
that this offspring would recombine with the remaining 



2A particles and undergo another local sequence of re- 
actions. And if in the meantime the 3^4 in a row have 
either branched or undergone a 3^4 — > 2A reaction, the re- 
maining particles would not have the time to go very far. 
Since everything takes place essentially locally, the mean- 
field rate equation should be valid [l8j yielding therefore 
a =1/3. 

To check this scenario I run simulations for the 3A — > 
4A, 3A — ► model, in which no such local sequence of 
reactions occur, since following the triplet annihilation 
the remaining single A can diffuse away at most. 

The dynamical rules are very similar to those of the 
3A — > 4A, 3A — > model. An elementary MCS consists 
of the following processes: 

(a) <-> 0A with probability D, 

(b) 3 A -> with probability p(l — D), 

(c) 3A -> 4A with probability (1 -p)(l - D), 

such that the reactions were allowed on the left or right 
side of the selected particle strings. Now the system size 
was 4 x L = 10 5 with periodic boundary conditions and 
the density of particles, singlets and triplets is followed 
up to 2 x 10 9 MCS from random initial state. As Figure 
01 shows the time evolution of pit) can be split into three 
parts (a) an initial slow regime (i<3x 10 3 MCS), (b) 
an intermadiate faster then mean-field regime (3 x 10 3 < 
t < 3 x 10 6 MCS), (c) a mean-field regime for t > 3 x 10 6 
MCS. One can see a level-off in the time evolution of 
pi 1 / 3 for p = p c = 0.11887(2). For this critical p value 
the density of triplets behaves in the same way in the long 
time limit (t >~ 3 x 10 6 MCS) corroboting the 4 = 1 
result of the previous section. 

Finally the static exponent (3 was determined in the 
active phase in the neighbourhood of the critical point 
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(p c = 0.11887). As Figure [3] shows the local slopes 
converge to the mean-field value again. By fit- 
ting with the form J2J one gets the following estimates: 
(3 = 0.99(5), x = 0.6(1). 



IV. CONCLUSIONS 

Large scale simulations for two different triplet mod- 
els: 3A -> 4A, 3A -> 2A and 3A -> 4A, 3A -> result in 
mean-field type of static and dynamic scaling behavior in 
one dimension. The simulations up to t = 2 x 10 9 MCS 
do not support the scenario according to which fluctua- 
tions are supressed and reactions take place locally hin- 
dering to see PCPD type of critical behavior. Since the 
triplet density dacays in the same way as the total den- 
sity, which is typical at the upper critical dimension the 
d c = 1 is concluded. Furthermore logarithmic correction 
to scaling is shown in case of the static order parameter 
exponent. On the other hand one can't see log. correc- 
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tions in the density decay, which may mean that these 
corrections are small or perhaps a next to the leading 
order correction term hinders to see it. Such correction 
term was found in a very recent field theoretical analysis 
of the 3 A — > (0, A, 2 A) models |l9j, which correspond to 
the dominant behavior in the inactive phase of the triplet 
model I studied. 

The contradiction with the results of the Langevin 
equation analysis of a bosonic triplet system is unre- 
solved. This may also mean that the two-species coupled 
model and the one-species model behave differently. 
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